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Abstract 

Recent developments of high-order CCM have been to extend existing formalism and codes to 
s > \ for both the ground and excited states, and independently to "generalised" expectation 
values for a wide range of one- and two-body spin operators. An advantage of the CCM is that the 
Goldstone linked-cluster theorem is obeyed at all levels of approximation and so it provides results 
in the infinite lattice limit N — > oo from the outset. However, recent results have also shown that 
the CCM can provide exact (symmetry-breaking) results for the spin-half linear-chain J1-J2 at the 
Majumdar-Ghosh point J2/J1 = 0.5 by identifying special solutions of the CCM equations for the 
usual Neel model state. Interestingly, the CCM provides exact (non-symmetry-breaking) results 
for systems in which small magnetic clusters become de-coupled from each other when the bonds 
connecting them tend to zero. These exact results involve the identification of "special solutions" 
of the CCM equations for the Neel state. An example of this is given by a spin-half system with 
nearest-neighbour bonds for an underlying lattice corresponding to the magnetic material CaV^g 
(CAVO) in which one of the two different types of bonds on the lattice tend to zero. Larger 
finite-sized systems may be considered by appropriate choice of the unit cell and the bonds on it. 
We show here that exact diagonalisation results for ground-state energy and excitation energy gap 
for the spin-half and spin-one linear Heisenberg model on chains of length up to N = 12 sites for 
s = 1/2 and N = 6 sites for s = 1 with periodic boundary conditions are reproduced exactly using 
high-order CCM via this "brute-force" approach; i.e., one in which none of the translational or 
point-group symmetries of the finite lattice are used. 
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I. INTRODUCTION 



The coupled cluster method (CCM) l|-[9| is a well-known method of quantum many- 
body theory (QMBT). The CCM has been applied with much success in order to study 
quantum magnetic systems at zero temperature (see Refs. |10h52||). In particular, the use 



of computer-algebraic implementations 



22 



26 



32 



52[ of the CCM has been found to be 



very effective with respect to these spin-lattice problems. Recent developments of high- 
order CCM formalism and codes have been to treat systems with spin quantum number of 
s > g for both the ground and excited states 52J. Furthermore, the ground-state formalism 



and codes may also be used directly to find "generalised" expectation values I 



These 



expectation values are defined for a wide range of one- or two-body spin operator that prior 
to the CCM calculation. 

Here we show how the consideration of previous results for exact results for the 



(symmetry-breaking) ID J1-J2 model at J2/J1 = 0.5 [50|,l54| and a (non-symmetry-breaking) 
nearest-neighbour CAVO model jso]] in the limits that various nearest-neighbour bond either 
go to zero or infinity leads on naturally to the treatment of finite lattices via the CCCM 
code 55] . This is achieved by the simple expedient of choosing the finite-lattice to be the 
fundamental unit cell and so this is a "brute-force" solution of the finite-lattice problem via 
high-order CCM. 



II. METHOD 



The details of the practical application of high-order coupled c 



uster method (CCM) 



formalism to lattice quantum spin systems are given in Refs. [22|, [26|, |32|, |38|, [52| and also in 
the appendices to this article. However, we point out now that the ket and bra ground-state 
energy eigenvectors, |\&) and of a general many-body system described by a Hamiltonian 
H, are given by 

H\*) = E g \V); {V\H = E g (V\. (1) 

Furthermore, the ket and bra states are parametrised within the single-reference CCM as 
follows: 
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(§| = ($|^e- 5 ; S=l + J2SiCT- (2) 

One of the most important features of the CCM is that one uses a single model or reference 
state |$) that is normalised. We note that the parametrisation of the ground state has the 
normalisation condition for the ground-state bra and ket wave functions = ($|$) = 

1). The model state is required to have the property of being a cyclic vector with respect to 
two well-defined Abelian subalgebras of multi-configurational creation operators {Cf } and 
their Hermitian-adjoint destruction counterparts {Cj = {Cfy}. The interested reader is 



52| for more information regarding how the CCM 



referred to the Appendices and to Ref. 
problem is solved for. 

Here, we use the Neel state as the model state for the antiferromagnetic Heisenberg model 
given by 

N 

H = JY. Bi-sj • (3) 

For the bipartite lattices, we perform a rotation of the local axes of the up-pointing spins 
by 180° about the y-axis. The transformation is described by, 

s x -> -s x , s y -> s y , s z -)• -s z . (4) 

The model state now appears mathematically to consist of purely down-pointing spins. In 
terms of the spin raising and lowering operators = s% ± is v k the Hamiltonian may be 
written in these local 



N 



h = -jj: 

(hi) 



s t4 + sTsJ + 2<s* 



(5) 



where the sum on again counts all nearest-neighbour pairs once on the lattice. 

The CCM formalism is exact in the limit of inclusion of all possible multi-spin cluster cor- 
relations within S and S, although this is usually impossible to achieve practically. Hence, 
we generally make approximations in both S and S. The three most commonly employed 
approximation schemes previously utilised have been: (1) the SUBn scheme, in which all 
correlations involving only n or fewer spins are retained, but no further restriction is made 
concerning their spatial separation on the lattice; (2) the SUBn-m sub-approximation, in 
which all SUBn correlations spanning a range of no more than m adjacent lattice sites are 
retained; and (3) the localised LSUBm scheme, in which all multi-spin correlations over all 
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distinct locales on the lattice denned by m or fewer contiguous sites are retained. Another 
important feature of the method is that the bra and ket states are not always explicitly 
constrained to be Hermitian conjugates when we make such approximations, although the 
important Hermann- Feynman theorem is always preserved. We remark that the CCM pro- 
vides results in the infinite-lattice limit iV — > oo from the outset. 

Key to understanding the application of the CCM to spin problems is the concept of the 
unit cell and the Bravais lattice. The unit cell contains a number of sites at specific positions 
(given by the "primitive" lattice vectors) that are replicated at all possible multiples of the 
Bravais lattice vectors. Thus, for example, we have a single site in the unit cell for the 
linear chain, say, at position (0,0) and a single Bravais lattice vector a = (1,0) T . The 
lattice is formed by translating the single site in the unit cell by all integer multiples of a. 
Two-dimensional lattices have two Bravais lattice vectors. For example the square lattice 
has a single site in the unit cell and the lattice vectors are a = (1,0) and b = (0, 1) T . 
The triangular lattice is given by vectors a = (1,0) T and b = (1/2, \/3/2) T and so on 
for other lattices. We see also that the basic building blocks of unit cell, Bravais lattice, 
and bonds/interactions in the Hamiltonian placed on the lattice gives us a broad canvas to 
work with. For example, we may form models that interpolate between different lattices 
(and even different spatial dimensions) by varying the strengths of various bonds that have 
been carefully placed with respect to the underlying lattice. Hence, the number of possible 
such quantum spin systems is enormous. Furthermore, the development in the number and 
complexity of these theoretical models is often driven by the magnetic materials studied in 
experiment. 

III. RESULTS 

A. The Spin-Half J\—Ji Model on the Linear Chain 

The Hamiltonian for this spin-half model has nearest-neighbour bonds of strength J\ and 
next-nearest-neighbour bonds of strength J 2 . We use a Neel model state in which nearest- 
neighbour spins on the linear chain are anti-parallel. We rotate the spin coordinates of the 
'up' spins so that notationally they become 'down' spins in these locally defined axes. The 
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relevant Hamiltonian in rotated coordinates is then given by 
H = 



(ij) «*,*» 



(6) 



where (i, j) runs over all nearest-neighbour sites on the lattice counting each pair once and 
once only and ((i, k)) runs over all next-nearest-neighbor sites on the lattice, again counting 
each pair once and once only. Henceforth we put J\ = 1 and consider J 2 > 0. 

The ground-state properties of this system have been studied using methods such as 



exact diagonalisations 



57 



62 



approaches 



r n 



DMRG 



17 



58 



60j, CCM 



421 ]. and field-theoretical 



60[ (see Refs. |60|,|61| for a general review). We shall not go into detail about 



this model here except to note that there are two degenerate simple exact dimer-singlet 
product ground states at the Majumdar-Ghosh point J2/J1 = 0.5. 

We now consider how this model can be treated at the Majumdar-Ghosh point J2/J1 = 
0.5by the CCM via the identification of a special dimerised solution of the CCM equations 
for a Neel model state. We use a doubled unit cell including two neighbouring sites for a 
spin-half system on the linear chain at points (0,0,0) and (1,0,0) and a single Bravais vector 
(2,0,0) T to take into account the symmetry breaking. There are thus two distinct types of 
two-spin nearest-neighbour ket-state correlation coefficients and again these are denoted as 
1S2 and 1S2 f° r ^2 = £2 Si a s t a s t a +i + ^2 Si b s t b s t b +i with respect to the CCM ket state of Eq. 
([2]), and where i a runs over all sites with odd-numbered indices and % runs over all sites with 
even-numbered indices. The exact ground state at J2/J1 = 0.5 is obtained by setting £2=1 
and all other coefficients equal to zero. This result is found to hold for all levels of LSUBm 
approximation. The exact ground-state energy of E g /N = — 0.375Ji is thus obtained at 
the point J2/J1 = 0.5, as expected. Furthermore, the sublattice magnetisation is found to 



be zero at this point at all levels of approximation. We find values [54J for the excitation 
energy gap of 0.35250, 0.34170, 0.30548, 0.28732, 0.27559, and 0.26760 at the LSUBm levels 
of approximation with m = {4,6,8, 10, 12, 14}. A simple extrapolation of these results in 
the limit m — > 00 using a quadratic function gives a value for the gap of 0.2310. This result 
is in agreement with results of exact diagonalisations that predict a gap of 0.234 63]. 



5 



FIG. 1: The CAVO lattice. The nearest-neighbour bonds that connect two sites on a four-site 
plaquette are shown by the solid lines and have a bond strength given by J (=1). The nearest- 
neighbour bonds that connect two sites on different plaquettes (dimer bonds) are shown by the 
dotted lines and have a bond strength given by J' . The unit cell of the lattice is shown by the 
square with the grey dashed lines. 

B. The J— J' Heisenberg antiferromagnet on the CAVO lattice 

In the previous section, we saw that the exact dimerised state at the Majumdar- Ghosh 
point may be obtained for the Neel model state by appropriate choice of ket-state correla- 
tion coefficients. However, this is a case in which the translational symmetry of the lattice 
(the linear chain) is spontaneously broken. We now wish to a case in which dimerised 
and plaquette exact solutions follow the symmetry of the Hamiltonian and so do not form 
symmetry-broken states. We consider an antiferromagnetic Heisenberg model in which the 
basic geometric unit cell contains four neighbouring lattice sites on the underlying crystal- 
lographic lattice of the magnetic material CaX^Og (CAVO), shown in Fig. [TJ There are 
two non-equivalent antiferromagnetic nearest-neighbour bonds J and J' belonging to dimers 
(J 7 ) and to four-spin plaquettes (J) respectively. The ground state of the quantum model 
depends on the ratio J' / J of the competing bonds. Using a unit cell as denned in Fig. [TJ 
the plaquette bonds J are inside the four-site unit cell and the dimer bonds J' connect sites 
in different unit cells. We note that this model is not frustrated but the two non-equivalent 
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nearest-neighbour bonds lead to a competition in the quantum system. Henceforth, we 
choose an energy scale such that J = 1. 

The four-site plaquettes in the unit cell become de-coupled only in the limit J'/ J — > 0. 
The ground state is a product of such four-site plaquette singlets in this limit. In the 
limit that J' /J — > oo dimers are formed on the J' bonds. To model such states using 
the CCM we start again from the Neel model state; namely, a state in which the spins 
on nearest-neighbour sites are anti-parallel. To create an exact plaquette-singlet product 
VBC ground state at J' /J = using the CCM we have to adjust the nearest-neighbour 
correlation coefficients «Sf and S\ and a single four-body plaquette correlation coefficient <Sf 
containing all four sites properly. (Note that S% represents those ket-state coefficients for the 
nearest-neighbour two-body cluster connecting sites on a plaquette indicated by the solid 
lines in Fig. [TJ whereas S\ represents those ket-state coefficients for the nearest-neighbour 
two-body cluster connecting sites on a dimer indicated by the dotted lines in the same 
figure. The coefficient S% represents those ket-state coefficients for the four-body cluster 
corresponding to a plaquette indicated by the solid lines in Fig. [TJ) Indeed, it is easy to 
show that setting the ket-state correlation coefficients S% and «S| to a value of 0.5 and all 
other ket-state correlation coefficients (including <S|) to zero the plaquette-singlet product 
valence-bond crystal state is obtained exactly. Furthermore, we are also able to reproduce 
exactly the dimer-singlet product ground state in the limit J' /J — > oo. In this limit, the 
nearest-neighbour ket-state correlation coefficient S\ on the dimer bonds (dotted lines in 
Fig. [TJ) has a value of one and all other coefficients (e.g., «Sf and Sf) are zero. 

An important point is that in the limits J' /J — > and \J'/J\ — > oo the system is 
comprised of independent clusters. However, the system is two-dimensional for all other 
values of J' /J. This system therefore "interpolates" between a zero-dimensional and two- 
dimensional lattice with the bond strengths J' /J. This model (and similar models) may 
therefore be used to study the differences between zero-dimensional and two-dimensional 
systems. They may also be used to investigate the effects of "linking" magnetic clusters 
in order to form an extended two-dimensional material. Such models are of interest in the 
subject of quantum computational. We remark that the CCM solution at intermediate values 
of J'/ J may be found by "tracking" the exact CCM solution by making small incremental 
changes in J'/ J from either of the exact limits J' / J — > and | J'/ J\ — > oo (see Ref. 50] for 
details). 
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FIG. 2: CCM results for the ground-state energy of the J~J' Heisenberg antiferromagnet on the 
CAVO lattice (with J = 1). 

We remark that the correct ground-state energies (with J = 1) are reproduced in the 
limits J'/ J and J'/ J — > oo, namely, E g /N = —0.5 and E g /N = —0.375 J', respectively. 
These results for the ground-state energy are shown in Fig. [21 and we see that LSUB8 and 
LSUB10 results agree extremely well over the region considered. Hence, we conclude that 
the CCM also provides excellent results for intermediate values of J' /J that interpolates 
between the two (exact) limits, J'/ J — > and J' / J — )■ oo. 

C. The Spin-Half and Spin-One Heisenberg Model on Finite-Sized Chains 

We have seen in the previous sections that an appropriate choice of the unit cell and values 
for the CCM correlation coefficients can reproduce exact results for Neel model states. This 
was observed, namely, for the symmetry-breaking solution to the J1-J2 model at J2/J1 = 0.5 
or for the CAVO model with two types of bonds in the limits J'/ J — > and \J'/J\ — > 00. 
We now consider if exact solutions for larger unit cells may be found by "brute force", i.e., 
without using any of the translational or point-group symmetries of the finite lattice. 

We start by defining the unit cell to be of size = {2,4,6,8,10,12}. We use the 
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Hamiltonian of Eq. fl3]) and we rotate the "up" spins as usual. The Hamiltonian in the 
new spin coordinates is then given by Eq. (jSJ). Hence, we form bonds (of strength J = 1) 
between all nearest-neighbour sites in the unit cell. However, we do not "link" the bonds 
by forming an intermediate bond between the unit cells, and so, in effect, the unit cells 
become isolated clusters. A technical point is that we set the single Bravais lattice vector to 
be of magnitude +N for this "brute force" approach; again, even though the clusters/unit 
cells are not linked by intermediate bonds. Furthermore, we may also provide an additional 
bond between the first and last sites (again of strength J = 1). In effect, this bond acts 
as a "boundary condition" that creates a "ring" of sites for the finite lattice. These results 
should be compared to those of ED with periodic boundary conditions. 

High-order CCM SUBn-m results for the antiferromagnetic Heisenberg model (with pe- 
riodic conditions) for finite-sized chains of length N with SUBn-m with n = m = N for 
s = 1/2 and n/2 = m = N for s = 1 in Tables HI and [H] respectively. We remark that CCM 
results for the ground-state energy and excitation energy gap agree to at least six decimal 
olaces with exact diagonalisations obtained using the SpinPack code of Joerg Schulenburg 
641 ] . We may also study the manner in which LSUBm results behave with increasing m for 
a set value of chain length N. Table HTT1 presents LSUBm results for the spin-half Heisenberg 
chain of length N = 12. The ground-state energies decrease monotonically with increasing 
LSUBm level of approximation, although no simple extrapolation "rule" may be seen. By 
contrast, the LSUBm data for the excitation energy gap is only monotonically decreasing up 
to LSUB10. Indeed, the gap for LSUB10 lies lower than that of the (exact) LSUB12 result. 
We note again that the LSUBm approximation for the spin-half system and the SUB2m-m 
approximation for the spin-one system reproduce exact results for the chains of length TV 
when we set m — N. 

The CCM results in Tables HI to IIIII arise from "special solutions" of the CCM equations in 
exactly the same manner as for the nearest-neighbour CAVO problem in the limits J' / J — > 
and | J'/ J\ — > oo. However, it is clearly more complicated in this case because we are dealing 
with unit cells (i.e., finite-sized ID lattices here) of greater size than 4 sites. Interestingly 
though, the CCM ground-state equations were found to converge readily to this "special 
solution" and for a wide range of starting values at all levels of SUBn-m approximation 
attempted here. For example, setting the initial ket-state correlation coefficients to a small 
(positive) non-zero value was found to work quite well. Indeed, the ground-state solution 
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TABLE I: CCM LSUBm results for the ground-state energy and excitation energy gap of spin-half 
antiferromagnetic Heisenberg chains of ./V length with periodic boundary conditions. 



N= 


CCM Eg/N 


CCM Gap 


4 


-0.5 


1 


6 


-0.46712928 


0.684740 


8 


-0.45638668 


0.522676 


10 


-0.45154464 


0.423239 


12 


-0.44894924 


0.355848 



TABLE II: CCM SUBn-m results (with n = 2N and m = N) for the ground-state energy and 
excitation energy gap of spin-one antiferromagnetic Heisenberg chains of N length with periodic 
boundary conditions. 



N= 


E g /N 


Gap 


2 


-1.263853992 


1.4876903286 


4 


-1.29781459 


1.038539433 


6 


-1.436237197 


0.720627363 



was found to be very stable, as is generally also the case for extended systems (N — > oo) in 
which the model state is known to be a "good starting point" . 

These results prove the principle that the CCM may be used to study finite-sized lattices. 



IV. CONCLUSIONS 



Again, we remark that recent developments of high-order CCM have been to extend 
existing formalism and codes to s > ~ for both the ground and excited states [52J, and 
independently to "generalised" expectation values for a wide range of one- and two-body 
spin operators 53]. We note that the CCM is that the Goldstone linked-cluster theorem is 
obeyed at all levels of approximation and so it provides results in the infinite lattice limit 
N — > oo from the outset. In this article, we have shown the exact ground state for the 
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TABLE III: CCM LSUBm results for the ground-state energy and excitation energy gap of spin-half 
antiferromagnetic Heisenberg chains of length N = 12 with periodic boundary conditions. 



LSUBm m = 


E g /N 


Gap 


4 


-0.421064786 


0.589781571 


6 


-0.42329432 


0.440170513 


8 


-0.423975557 


0.375801129 


10 


-0.424415368 


0.348833726 


12 


-0.448949243 


0.355847514 



J\-Ji model at J2/J1 = 0.5 may be reproduced using a Neel model state. Using a similar 
model state in 2D, exact results for a nearest-neighbour CAVO model were reproduced using 
high-order CCM in the limits J'/ J — > (4-site plaquette) and J' /J — > 00 (2-site dimers). 
These results lead on naturally to a "brute-force" approach for solving finite-sized lattices, 
i.e., without using any of the translational or point-group symmetries of the finite lattice. 
Indeed, we have shown here that high-order CCM SUBn-m results for the antiferromagnetic 
Heisenberg model (with periodic conditions) for finite-sized chains of length A" with SUBn- 
m with n = m = A r fors = l/2 and n/2 = m = N for s = 1 were found to agree with 
exact-diagonalisation (ED) results for the ground-state energy and excitation energy gap 
agree to at least six decimal places. 

We note that results for the case of finite clusters for s = 1 were obtained via new high- 
order formalism for the CCM excited state outlined in Refs. 52j, |54j. Furthermore, new 
code 55[ has been written to implement this new formalism for the excited state for s > 1. 
The agreement between ED results for the ground-state energy and excitation energy gap 
for chains of up to A" = 6 for s = 1 is an excellent test of the validity of this new code. These 
solutions for the finite chains were found to be stable numerically. Hence, we have "proven 
the principle" that the high-order CCM code may be used directly to study (relatively small- 
sized) finite lattices by a somewhat "brute-force" approach. Indeed, this approach is still 
somewhat inefficient in comparison to exact diagonalisations (ED) because ED results use 
translational and point-group symmetries of the finite lattice in order to reduce the size of 
the matrix to be diagonalised, and so ED may to go to much larger lattice sizes. 
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For periodic boundary conditions, the manner in which this is achieved for ED is by 
identifying common states for a given fc-value via a complex phase factor e _lk r for states 
that a related by a translational vector r on the finite lattice (again: note periodic boundary 
conditions are assumed). The Hamiltonian then links only those states of common k. Point- 
group (PG) symmetries/permutations of indices may then used to form states with real 
components only, thus simplifying the computational problem. In principle, we ought to 
be able to employ the finite lattice translational and point-group symmetries analogously 
in order to simplify the finite-size problem with periodic boundary conditions also for the 
CCM. However, it is unclear how one might do this in practice for the exponentiated S in 
the ground ket and bra states. The use of translational and PG symmetries for finite-lattice 
CCM will be the subject of future research. 

Despite the fact that translational and PG symmetries for the finite lattice were not 
used directly here, we were still able to treat finite-lattices of size N = 12 for the spin- 
half case using high-order CCM with only relatively meagre computational resources. (A 
MacBook with a "Core Duo" processor and 1 GIG RAM was used in this case was used 
to carry out this calculation.) Still larger lattices are possible using the CCCM code 55], 
which has been implemented to work in parallel on a cluster of processors. We note that 
the addition of intermediate bonds that link the isolated clusters in order to form extended 
lattices of infinite numbers of sites requires only relatively small increases in computational 
effort for high-order CCM compared to treating the case of isolated magnetic clusters alone. 
Excellent results were seen here for such a model that interpolated between finite clusters 
and an extended lattice, namely, for the n.n. CAVO model (e.g., see the results for the 
ground state energy in Fig. [2]). High-order CCM might provide a good choice for the study 
of a whole range of such "interpolating" models between finite clusters and infinite lattices. 
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